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Abstract. A study of the O/Fe ratio in metal-poor main sequence and subgiant stars is presented using the 
[0 i] 6300 A line, the Ol 7774 A triplet, and a selection of weak Fell lines observed on high-resolution spectra 
acquired with the VLT UVES spectrograph. The [O i] line is detected in the spectra of 18 stars with —0.5 < 
[Fe/H] — 2.4, and the triplet is observed for 15 stars with [Fe/H] ranging from —1.0 to —2.7. The abundance 
analysis was made first using standard model atmospheres taking into account non-LTE effects on the triplet: 
the [Ol] line and the triplet give consistent results with [O/Fe] increasing quasi- linearly with decreasing [Fe/H] 
reaching [O/Fe] ~ +0.7 at [Fe/H] = —2.5. This trend is in reasonable agreement with other results for [O/Fe] 
in metal-poor dwarfs obtained using standard atmospheres and both ultraviolet and infrared OH lines. There is 
also broad agreement with published results for [O/Fe] for giants obtained using standard model atmospheres 
and the [O i] line, and the OH infrared lines, but the O I lines give higher [O/Fe] values which may, however, fall 
into place when non-LTE effects are considered. When hydrodynamical model atmospheres representing stellar 
granulation in dwarf and subgiant stars replace standard models, the [O/Fe] from the [Ol] and Fell lines is 
decreased by an amount which increases with decreasing [Fe/H]. These 3D effects on [O/Fe] is compounded by 
the opposite behaviour of the [Ol] (continuous opacity effect) and Fell lines (excitation effect). The [O/Fe] vs 
[Fe/H] relation remains quasi-linear extending to [O/Fe] ~ +0.5 at [Fe/H] = —2.5, but with a tendency of a 
plateau with [O/Fe] ~ +0.3 for -2.0 < [Fe/H] < -1.0, and a hint of cosmic scatter in [O/Fe] at [Fe/H] ~ -1.0. 
Use of the hydrodynamical models disturbs the broad agreement between the oxygen abundances from the [O i] , 
O I and OH lines, but 3D non-LTE effects may serve to erase these differences. The [O /Fe] values from the [O i] 
line and the hydrodynamical model atmospheres for dwarfs and subgiant stars are lower than the values for giants 
using standard model atmospheres and the [O i] and O I lines. 
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1. Introduction 

The oxygen abundance of metal-poor stars is intimately 
linked to several outstanding questions in astrophysics. 
Questions of current interest include the following three: 
How old are the oldest (i.e., the most metal-poor) stars? 
Can the observed chemical evolution of oxygen in the early 
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Galaxy be interpreted using current predictions for the 
yields from metal-poor massive stars exploding as super- 
novae? What conclusions may be drawn from the abun- 
dance ratios of lithium, beryllium, and boron relative to 
oxygen about the synthesis of these three light elements 
by collisions involving high energy particles in interstellar 
circumstellar gas? For these and other reasons, accurate 
determinations of the growth of the oxygen abundance of 
the Galaxy are desirable. 

Oxygen, the most abundant element after hydrogen 
and helium, affects the opacity and, hence, the structure 
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and evolution of low mass stars. Globular clusters are con- 
sidered to contain some of the oldest stars in the Galaxy. 
Oxygen affects the shape of a cluster color-magnitude di- 
agram, and, hence, a cluster's age as it is determined 
from the fit of a theoretical to the observed diagram. 
VandenBerg & Bell (2001) estimate that 'each increase of 
0.3 dex in [O/Fe] implies a 1 Gyr age reduction at a fixed 
turnoff luminosity'. In the case of the most metal-poor 
stars, current estimates of the [O/Fe] differ by as much as 
0.7 dex or ages differ by 2 to 3 Gyr, a not insignificant 
difference for cosmologists. 

Compositions of metal-poor unevolved stars are pre- 
sumed to reflect the composition of the gas from which 
they formed. In the case of the most metal-poor stars, 
that gas may have been primordial gas mixed with ejecta 
from just one or two Type II supernovae. Hence, the pre- 
dicted compositions of the ejecta may be tested against 
observed compositions. Predictions appeared not to have 
a problem in matching observations that suggested [O /Fc] 
~ 0.4, as well as [a/Fe] ~ 0.4, where a denotes Mg, Si, 
and Ca. With some observers now reporting much higher 
oxygen abundances, [O/Fe] ~ 1 for very metal-poor stars, 
serious challenges may face the modelers of the Type II 
supernova (see particularly Goswami & Prantzos 20001 ). 
Recent claims that hypernovae, i.e. supernovae with very 
high explosion energies, at low metallicities could account 
for the high values of [O/Fe] (Israelian et al. 2001) ap- 
pear not to be supported by the most recent theoretical 
investigations of such events. On the contrary, detailed nu- 
cleosynthetic calculations reveal significantly lower O/Fe 
ratios than in normal core-collapse supernovae (Nakamura 
et al. pOOl]). 



Beryllium is identified as solely a product of spalla- 
tion of nuclei in interstellar or circumstellar gas traversed 
by energetic particles: for example, the process p + O — > 
splits 9 Be nuclei from the oxygen 16 nucleus. Spallation 
also makes a contribution to the synthesis of lithium and 
boron. Under quite simple but reasonable assumptions 
concerning the nucleosynthesis of oxygen and the origin 
of the energetic particles (here referred to as cosmic rays, 
or CR for short), if beryllium synthesis is controlled by 
a CR proton (or a-particle) colliding with an interstellar 
oxygen nucleus, one expects the beryllium abundance to 
grow as the square of the oxygen abundance. If the roles 
are reversed and a cosmic ray oxygen nucleus collides with 
an interstellar proton (or a-particle) , the beryllium and 
oxygen abundances are linearly related. In the present in- 
terstellar medium, this latter process is a minor contrib- 
utor to beryllium synthesis because the fragments of the 
spallation of the oxygen nuclei have high energies and are 
themselves spallated as they are thermalized. In the for- 
mer process, the fragments have much lower energies and 
a higher survival rate in thermalization. These relative 
roles may be reversed in the early Galaxy when the in- 
terstellar oxygen abundance was much lower than in the 
interstellar medium at present. By determining the Be/O 
ratio for metal-poor (and metal-rich) stars, it is possible 
to gain insight into the spallation processes. Uncertainty 



about the oxygen abundance necessarily compromises this 
exercise. Substitution of iron for oxygen is unsatisfactory 
because very little beryllium is produced from iron, and 
also the predicted oxygen to iron ratio in ejecta of Type 
II supernovae is uncertain. 

Spectra of metal-poor stars offer several potential indi- 
cators of the oxygen abundance: ultraviolet OH lines, the 
[O i] 6300 A and 6363 A lines, the O I triplet at 7774 A, and 
infrared vibration-rotation lines from the OH molecule's 
ground state. A selection of these lines is available for both 
metal-poor dwarfs and giants. A problem of great current 
concern is that the various indicators may give different re- 
sults. An abbreviated discussion suffices here, and a fuller 
discussion is postponed to Section 5. 

The problem became strikingly evident with the pub- 
lication of analy ses of the ultraviolet O H line s in dwarfs 
(Israelian et al. 1998; Boesgaard et al. 1999). These in- 



dependent but quite similar analyses showed that [O/Fe] 
rose linearly with [Fe/H] from zero at [Fe/H] = to about 
1.2 at [Fe/H] = —3.0, about the metallicity of the most 
metal-poor stars studied. Accounting for non-LTE effects 
on the reference Fe I lines diminishes the slope somewhat 
to —0.33 but still with a monotonic trend (Israelian et 
al. 2001). The linear increase to lower metallicities antic- 
ipated by Abia & Rebolo (1989) from the Oi lines was 
confirmed by some but not all analyses of these lines. A 
quite different relation for [O/Fe] versus [Fe/H] had been 
derived prior to the OH analyses from the [Oi] lines in 
giants: [O/Fe] ~ 0.4 independent of [Fe/H] for giants with 
[Fe/H] < -1 (Barbuy |l98S| ). At [Fe/H] = -3, the OH in 
dwarfs and the [O i] in giants results for [O /Fe] differ by 
about 0.6 dex. 

A strong a priori claim may be advanced that the 
[O i]-based analysis is the least sensitive to systematic er- 
rors: the lines are formed in or very close to local ther- 
modynamic equilibrium (LTE); essentially all of the oxy- 
gen atoms in the dwarf and giant atmospheres are in the 
ground configuration providing the lower and the upper 
level of the forbidden lines. Even the sensitivity of the 
[O i]-based results to the assumed surface gravity is effec- 
tively mitigated by determining the iron abundance from 
Fen lines. This pairing also ensures the non-LTE effects 
affecting the iron abundance are minimal; stronger non- 
LTE effects are predicted for Fel lines (e.g. Thevenin & 
Idiart 1999). Dependence of [O/Fe] on effective temper- 
ature is also quite insignificant when the combination of 
[O i] and Fell lines is used, whereas the dependence on the 
atmospheric structure in general cannot be neglected, as 
will be discussed later, but at least it is smaller than in 
the case where OH lines are used. 

By sharp contrast, the number of OH molecules (and of 
O I atoms in the lower state of the 7774 A transitions) is a 
very small fraction of the total number of oxygen atoms, 
a fraction that is highly sensitive to the adopted values 
of effective temperature and surface gravity. Furthermore, 
LTE is far from gu arant eed for either the OH UV lines 
(Hinkle & Lambert |l975| ; Asplund & Garcia Perez |200l|) 
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or the O I 7774 A lines (Kiselman 2001 and references 
therein) . 

The view has been expressed that the [O i]-based re- 
sults for the giants are in error on two accounts (Takada- 
Hidai et al. 2001; Israelian et al. 2001): (i) oxygen is re- 
duced in the atmospheres of giants because the convective 
envelope has mixed very oxygen-poor material to the sur- 
face from the ON-cycled interior, and (ii) the [Oi] line is 
partially filled in by emission. These are suspicions with- 
out firm supporting observational evidence. Reduction of 
oxygen by mixing must lead to a severe enhancement of 
the nitrogen abundance, a result not suspected from exam- 
ination of spectra and not yet confirmed by quantitative 
analysis. If emission is the culprit, is it not surprising that 
the [O i] lines have never been seen in emission in metal- 
poor giants, and that the absorption line profile has not 
been reported as distorted or shifted in radial velocity? 

A new way to approach the problem of the oxygen-to- 
iron ratio in metal-poor stars is to determine the abun- 
dances from observations of the [Oi] and Fen lines in 
metal-poor ([Fc/H] < —1) main sequence and subgiant 
stars. Mixing in these stars is not predicted. Emission in 
the [O i] lines would seem highly improbable. The diffi- 
culty is that the [O i] lines are very weak in the spectra 
of such stars, and up to now only 4 dwarfs (Spiesman 
& Wallerstein |1 99 1| ; Spite & Spite |1991|) and 2 subgiants 
(Fulbright & Kraft [199$ Cayrel et al. 200l| ; Israelian et 
al. pOOl ) have been studied. It is our goal in this paper 



to provide a more extensive study of the [O i] 6300 A, the 
stronger of the 2 lines, in metal-poor main sequence and 
subgiant stars based on new VLT/UVES spectra. At the 
same time, we determine the iron abundance from weak 
Fen lines, which makes the derived O/Fe ratio practically 
independent of non-LTE effects, uncertainties in the as- 
sumed surface gravity, and only weakly dependent on er- 
rors in the effective temperature, T e ff. Furthermore, we 
compare the [O i] based abundances with non-LTE oxy- 
gen abundances derived from the Oi 7774 A triplet. 

2. Observations and data reduction 

The oxygen and iron abundances have been derived from 
spectra obtained with UVES, the Ultra violet and Visual 
Echelle Spectrograph (Dekker et al. [2000D at the ESO VLT 
8m Kueyen telescope, during three observing runs: April 
8-12, 2000, July 22-24, 2000 and April 10-12, 2001. The 
main purpose of the April observations was to determine 
beryllium abundances from the Ben lines at 3130 A as ob- 
served in the blue arm of UVES. Using a dichroic filter the 
spectral range 4800-6800 A region was observed simultane- 
ously in the red arm of UVES with a resolution of 55 000 
and about 4 pixels per spectral resolution element. The 
July 2000 observations were primarily conducted to de- 
termine the 6 Li/ 7 Li isotopic ratio from the profile of the 
Li I resonance line at 6708 A. The spectra were obtained 
with a narrow slit (0.3 arc sec) in combination with an 
image sheer and have a resolution of 120000 with two 
pixels per spectral resolution element. Only the red arm 



of UVES was used and the spectral ranges covered were 
6000 - 7000 A with the EEV CCD and 7100-8200 A with 
the MIT chip. 

The spectra have been reduced using MIDAS for the 
April spectra and IRAF for the July spectra. Extensive 
testing has been performed between these two data reduc- 
tion packages, showing that the final spectra are indeed of 
similar quality. In both cases the reduction includes defi- 
nition of the echelle orders, background subtraction, flat- 
field division, extraction of orders, wavelength calibration 
and continuum normalization. The S/N of the resulting 
spectra are very high: 300-400 per pixel for the April spec- 
tra and 400-600 for the July spectra. The continuum was 
fitted using cubic spline functions with a wavelength scale 
of about 5 A. As the spectrum of a metal-poor star con- 
tains many good continuum regions this ensures that the 
continuum is well defined except in the vicinity of broad 
spectral features like the Ha line and the bandhead of 
telluric O2 lines. 

Depending on the radial velocity of the star, the [0 1] 
line may be blended by telluric O2 and H2O lines. In order 
to be able to remove these lines, several rapidly rotating 
B-type stars were observed. The IRAF task telluric was 
then used to optimize the fit between the telluric lines 
of the program star and the B-type star by allowing a 
relative scaling in airmass and a wavelength shift. The 
procedure worked well for the July 2000 observations due 
to the high resolution and a stable spectral PSF of UVES 
as defined by the output from the image sheer and the 0.3 
arcsec entrance slit, but was less succesfull for the April 
observations, apparently because the PSF in this case is 
defined by a convolution of the entrance slit (0.7 arcsec) 
and the variable seeing profile of the star. Hence, we were 
not able to use the [0 1] line in the April spectra if it was 
significantly overlapped by telluric lines. 

Figure yj shows some representative spectra. The first 
two are from the April 2000 observations and have S/N 
close to 400, whereas the last two are from July 2000 and 
have S/N ~ 600. We note that in the case of HD 140283 
the [Oi] line was blended by a weak telluric H2O line 
(W ~ 1mA). Care was taken to remove this line by a 
proper scaling of the B-type calibration spectrum based 
on numerous other H2O lines, but the presence of the 
blending H2O line adds to the uncertainty of the oxygen 
abundance derived for HD 140283. 

3. Stellar parameters 

On the basis of available Stromgren photometry and other 
literature data the programme stars were selected to lie 
rather close to the turnoff region of halo stars and to span 
the metallicity range —3.0 < [Fe/H] < —0.6. The oxygen 
abundance could not be determined for all stars because 
either the [0 1] 6300 A line was too weak to be detected or 
it was blended with telluric lines. Still, we want to derive 
stellar parameters, T c g, logg and [Fe/H] for all program 
stars. These parameters will also be needed in connection 
with the study of Li and Be abundances. 
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Fig. 1. Representative spectra in the 6300 A region ar- 
ranged in order of decreasing strength of the [Oi] line. 
The spectra are corrected for Doppler shifts to the labo- 
ratory wavelengths of the stellar lines. In the spectrum of 
HD 111980 a telluric O2 line is blending the Sil line and 
in the case of HD 3567, the atmospheric emission line of 
[Oi] is seen to the right of the Fel line. 



3.1. Effective temperature 

T e ff was determined from the b — y and V — K colour in- 



dices using the IRFM calibrations of Alonso et al. (1996a). 
Stromgren uvby-/3 photometry was taken from Schuster & 
Nissen ( |1988| ) for the large majority of the stars supple- 
mented with unpublished photometry of Schuster et al. 
( pOO^ ) for the remaining stars. The (b — y)o calibration of 
Schuster & Nissen ( 1989] ) ( includ ing a zero-point correc- 
tion of +0.005 mag, Nissen 1994 ) was used to derive the 
interstellar reddening excess. If E(b — y)> 0.015, the red- 
dening is considered significant and the V magnitude as 
well as the mi and the c\ indices are corrected according to 



the relations V = V - A.3E(b - y) , m = m x + 0.3E(b - y) 
and Co = c% — 0.2E(b — y). These values and the reddening 
estimates are given in Table 1. 



The K photometry was taken from Carney (1983a), 
Alonso et al. ( [1994) ) and The Two M icron All Sky Survey 
(2MASS) data base (Finlator et al. |2000[ ). Furthermore, 
Dr. A. Alonso kindly supplied unpublished K magnitudes 
for a few additional stars and Dr. G. Simon sent us K pho- 
tometry for a number of southern stars observed incon- 
nection with the DENIS project (Epchtein et al. 1999 ). 
In general, the photometry of the various sources agree 
within a typical scatter of ±0.04 mag for stars in common 
and no obvious systematic differences are seen. A straight 
mean value of K from the various sources has therefore 
been adopted. Table 1 lists the corresponding V — K in- 
dex corrected for interstellar reddening according to the 
relation E (V - K) = 2.7 E(B -V) = 3.8E(b - y) (Savage 
& Mathis P.979D if E(b - y) > 0.015. 

When deriving T e g from b — y it was noted that the cal- 
ibration equation of Alonso et al. ( 1996a , Eq. 9) contains a 
quadratic term in [Fe/H], which introduces a tendency for 
T e ff to diverge towards high values at very low metallici- 
ties. As pointed out by Ryan et al. (1999) this behaviour 
(see their Fig. 5d) is likely to be a numerical artifact of 
the calibration. We have investigated the problem by com- 
paring T c ff derived from the Alonso et al. calibration with 
original IRFM temperatures from Alonso et al. (1996b) 
for a sample of 49 stars with Stromgren photometry and 
[Fe/H] < —1.5. Indeed, there is a trend that stars with 
[Fe/H] < —2.0 are assigned too high a T c s from the Alonso 
et al. calibration. To avoid this problem we have simply 
adopted a lower limit of [Fe/H] = —2.1 in their Eq. (9). 
The T cB calibration ofV — K (Alonso et al. [19964 E q- 7 ) 
has also a square term in [Fe/H]. The coefficient is, how- 
ever, very small and the term is quite insignificant even at 
the lowest metallicities, [Fe/H] ~ —3.0. 

The typical observational error of b — y is 0.007 mag, 
which corresponds to an error of ±50 K in T g. The error 
of V — K may be as high as 0.05 mag, but due to the 
high sensitivity of V — K to T G g this also leads to an 
error of ±50 K only. An independent estimate of the error 
may be obtained by comparing T c g derived from b — y 
and V — K, respectively Table 1 lists the data and Fig. 
U shows the comparison. As seen, T e g(b — y) is on the 
average 60 K higher than T c s(V — K). We have not found 
any explanation of this offset, but we note that there are 
no significant correlations between the deviation and T e ff, 
absolute magnitude, or [Fe/H]. Hence, we adopt the mean 
value of T c ff(& — y) and T e g(V — K). For the six stars 
missing V— K, a value T e g (b — y) — 30 K has been adopted. 

The rms scatter of the T e g difference shown in Fig. 
is 68 K. This corresponds to a statistical error of about 
±50 K on T e g (b — y) and T e g(V — K) assuming that they 
have the same error. If both b — y and V — K are avail- 
able we should then be able to reach a statistical error 
of T e ff of about ±35 K, which is co mparable to the ac- 
curacy claimed by Ryan et al. ( 1999 ) in their study of Li 
abundances of metal poor halo stars. This is, however, an 
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Table 1. Stromgren photometry, colour excess, V — K index, T e s from b — y and V — K, spectroscopic value of [Fe/H] , 
and absolute magnitudes derived from the Stromgren photometry and from the Hipparcos parallax including an error 
corresponding to the parallax error. If E(b — y) > 0.015, the V magnitudes and the photometric indices have been 
corrected for interstellar absorption. For the two SB2 stars the photometric value of [Fe/H] is given. 
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337 


2 


597 





007 






6144 




-1.94 


3 


55 


3 


43 


o 


59 


HD 298986 


10 


049 





324 





080 





301 










1 


312 


6081 


6061 


-1.35 


4 


35 


4 


48 





40 


HD 338529 


9 


370 





308 





045 





366 


2 


600 





006 


1 


260 


6278 


6200 


-2.34 


3 


22 


3 


57 





46 


CD-30° 18140 


9 


859 





302 





053 





340 


2 


606 





021 


1 


214 


6288 


6260 


-1.91 


3 


93 


4 


18 





46 


CD-35° 14849 


10 


568 





321 





040 





293 


2 


603 





010 






6147 




-2.34 


4 


24 










CD-57°1633 


9 


527 





343 





113 





271 










1 


418 


5944 


5873 


-0.89 


4 


66 


4 


67 





19 


G 013-009 


9 


916 





292 





054 





369 


2 


609 





019 


1 


219 


6402 


6275 


-2.30 


3 


68 


3 


71 





59 


G 020-024 


10 


764 





278 





059 





357 


2 


627 





084 


1 


132 


6469 


6411 


-1.78 


4 


33 










G 075-031 


10 


521 





333 





093 





330 


2 


597 


-0 


005 


1 


421 


6035 


5870 


-1.04 


3 


96 










G 088-032 


10 


780 





311 





039 





356 


2 


591 


-0 


002 


1 


235 


6251 


6252 


-2.41 


3 


26 










G 126-062 


9 


478 





330 





063 





327 


2 


588 


-0 


004 


1 


455 


6058 


5824 


-1.58 


3 


75 


4 


11 





37 


G 183-011 


9 


861 





319 





050 





310 


2 


602 





006 


1 


354 


6170 


6015 


-2.10 


4 


08 










G 271-162 


10 


264 





306 





055 





355 


2 


602 





020 


1 


275 


6289 


6170 


-2.31 


3 


61 










LP 815-43 


10 


774 





272 





054 





376 


2 


623 





032 


1 


140 


6558 


6471 


-2.72 


4 


17 











SB2 



SB2 



optimistic estimate of the error. Several of the stars, espe- 
cially the most metal-poor ones, are distant enough to be 
significantly reddened, and considering that the reddening 
estimate is hardly better than E(b — y) ~ 0.01 an error of 
70 K is introduced in both T cS (b - y) and T cS (V - K). 
Furthermore, T e g may be affected by systematic errors, 
which could depend on [Fe/H] and the position of the star 
in the HR diagram. 



3.2. Surface gravity 

The surface gravity was determined from the fundamental 
relation 



log- 



50 



log 



M 



M 







0A(M bol - M bolB ) 



4 log 



T Pl 



T 



ff,0 



(1) 



where M. is the mass of the star and M^oi the absolute 
bolometric magnitude. 

The absolute visual magnitude My was determined 
from a new cali bratio n of the Stromgren indices derived by 
Schuster et al. ( |2002 ) on the basis of Hipparcos paral laxes , 
and also directly from the Hipparcos parallax (ESA 1997 ) 
if available with an accuracy cr(7r)/7r < 0.3. Column 12 
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difference between effective temperatures de- 
y and V - K vs. T eB (V - K). 



o 
o 



it 



Q: 
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3.0 



4.5 
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Fig. 3. The difference between the absolute visual mag- 
nitude determined from Stromgren photometry and from 
the Hipparcos parallax. The error bars correspond to the 
parallax error and an estimated error of My, phot of 0.15 
mag. 



and 13 of Table 1 lists the photometric and the parallax 
based values of My and in Fig. [|the values are compared. 
As seen the agreement is quite satisfactory considering 
the estimated error bars. For two stars, HD 132475 and 
HD 189558, Mv, par has a 3a deviation from My.phot of 
the order —0.8 mag, which suggests that these stars are 
binaries consisting of two almost equal components. There 
is, however, no indication of double or asymmetric lines 
in their spectra, so we reject the binary hypothesis and 
consider the deviation in Fig. || as accidental. Two other 
stars, HD 99383 and HD 116064, have double lines, clearly 
indicating the presence of a fainter and cooler component 
to the main star. These spectroscopic binaries (SB2) stars 
were excluded from further analysis. 

In calculating logg we adopted the mean value of 
My.phot and My. par if the parallax is available, otherwise 
My, phot- The bolometric correction is taken from Alonso 



lating in the My- log T e ff diagram between the a-element 
enhanced evolutionary tracks of VandenBerg et al. ( 2000 ) . 
The internal precision of the mass is estimated to be better 
than ±0.05 Mq. 

The error of the derived value of logg is dominated 
by the error of M v ; a(M v ) = ±0.20 corresponds to 
er(logg) = ±0.08 dex. In addition there could be a small 
systematic error of log g but in view of the good overall 
agreement between Mv,phot an d My tPar this systematic 
error is probably less than ±0.10 dex. Altogether we esti- 
mate that the error of logg is on the order of ±0.15 dex. 

The final values of T c ff and log g are given in Table 2 
together with values for [Fe/H] and the microturbulence 
as derived from Fen lines (Sect. 4). In addition, the abso- 
lute magnitude and the mass are given. As the calibration 
of T e g and hence also logg depend somewhat on [Fe/H] 
the determination of the three parameters is an iterative 
process. 

4. O and Fe abundances 



The determination of abundances is based on a-element 
enhanced ([a/Fe]= ±0.4, a = C, O, Ne, Mg, Si, S, Ca, 
and Ti) ID model atmospheres with Tcff, logg and [Fe/H] 
values as given in Table 2 and a microturbulence of 1.0 
km s -1 . The models were computed with the MARC S code 
using updated continuous opacities (Asplund et al. 1997 ) 
and including UV line blanketing by millions of absorption 
lines. LTE is assumed both in constructing the models and 
in deriving abundances. We first discuss the determination 
of Fe abundances from Fe n lines and then O abundances 
derived from the [O i] 6300 A line and the O I triplet at 
7774 A. 



4.1. Iron abundance 

By inspecting the spectra of typical program stars as well 
as the solar spectru m, 13 a pparently unblended Fe n lines 
from Biemont et al. ( 1991 , Table 2) were selected as suit- 
able for determining the Fe abundance. The lines are listed 
in Table 3. The logg/ values given are those of Biemont et 
al., i.e. based on theoretical transition probabilities and a 
normalization (—0.05 dex) from lifetime measurements. 
The lifetime measurements reported by Biemont et al. 
agree well with the re-analysis by Allende Prieto et al. 
( 2002 ) of published experimental data, including more re- 
cent results. 

Using the MARCS solar model and assuming a micro- 
turbulence of 1.15 km s _1 the equivalent widths of the lines 
as measured in the solar flux spectrum of Kurucz et al. 



(1984) were first analyzed. The results are given in Table 3. 
The average logarithmic solar abundance, loge(Fe) = 7.53, 
is close to the commonly adopted meteoritic value of 7.50 



et al. ( 1995 Table 4). The mass was derived by interpo- 



(Grevesse & Sauval 1998), but the dispersion is quite high 
(0.08 dex). We note that there is no significant correlation 
between the Fe abundance and the equivalent width of 
the lines suggesting that the chosen value of the microtur- 
bulence is approximately correct. Furthermore, it is seen 
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Table 3. List of Fen lines used for determining the iron abundances of the stars. Equivalent widths as observed 
in the solar flux spectrum and in the spectra of HD 189558, HD 160617 and HD 140283 are given together with the 
derived abundances for each line. 



AVavelengt h 
A 


exc.pot. 
eV 


lot? (7 f 


Solar flux 
W(mk) loge(Fe) 


HD 189558 
W(mA) [Fe/H] 


HD 160617 
W(mk) [Fe/H] 


HD 140283 
W(mk) [Fe/H] 


4993.35 


2.81 


-3.54 


40.5 


7.33 


13.5 


-1.10 


3.8 


-1.81 






5100.66 


2.81 


-4.19 


19.6 


7.46 


4.5 


-1.12 


1.2 


-1.81 






5197.58 


3.23 


—2.28 


08.4 


1 CO 

/ .02 


trn o 

oo.o 


1 1 9 

— l.lo 


zy. / 


— 1.81 


1 A O 


O A *7 

— LA 1 


5234.62 


3.22 


-2.20 


90.5 


7.46 


56.4 


-1.09 


33.4 


-1.76 


12.7 


-2.39 


5325.55 


3.22 


-3.27 


43.0 


7.50 


15.1 


-1.09 


4.9 


-1.76 






5414.08 


3.22 


-3.80 


27.0 


7.65 


7.2 


-1.09 


2.5 


-1.68 






5425.26 


3.20 


-3.42 


44.0 


7.66 


14.6 


-1.14 


4.4 


-1.84 


1.4 


-2.42 


6084.10 


3.20 


-3.86 


21.3 


7.52 


5.5 


-1.07 


1.4 


-1.80 






6149.25 


3.89 


-2.77 


38.0 


7.51 


11.8 


-1.10 


3.9 


-1.76 


1.2 


-2.33 


6247.56 


3.89 


-2.38 


56.2 


7.53 


23.4 


-1.12 


8.7 


-1.80 


2.1 


-2.49 


6416.93 


3.89 


-2.79 


42.4 


7.62 


11.8 


-1.19 


4.2 


-1.83 






6432.68 


2.89 


-3.76 


43.0 


7.65 


14.8 


-1.13 


4.6 


-1.81 


1.3 


-2.46 


6456.39 


3.90 


-2.13 


66.5 


7.50 


33.0 


-1.09 


14.2 


-1.77 


4.3 


-2.39 


Average: 
Dispersion: 








7.53 
0.08 




-1.11 
0.03 




-1.79 
0.04 




-2.42 
0.06 



that when determining differential Fe abundances for each 
line for three representative stars, HD 140283, HD 160617 
and HD 189558, an amazingly small line-to-line scatter of 
[Fe/H] (0.03 - 0.06 dex) is obtained. This suggests that the 
scatter of the solar loge(Fe) values is mainly due to errors 
in the gf- values. By working differentially with respect to 
the Sun this problem is avoided. For the whole set of pro- 
gram stars the scatter in [Fe/H] ranges from 0.03 to about 
0.08 dex, and since at least four Fell lines are measured 
even in the most metal-poor stars (see Table Al [|), the 
statistical error of [Fe/H] is less than 0.05 dex. 

In calculating abundances from the Fell lines wc 
adopted the Unsold (1955) approximation to the Van der 
Waals interaction constant with an enhancement factor 
E 1 — 2.0. We note that a decrease of E 7 to 1.0 leads 
to an increase of the solar iron abundance of +0.04 dex, 
whereas the Fen lines in the stars are weak enough to 
make the derived iron abundance practically independent 
of Ey. Hence, the adopted value of E~ has no more than 
a small systematic effect on the metallicity scale of the 
metal-poor stars and the derived [O/Fe] values. 

In the more metal-poor stars, [Fe/H] < —1.5, the Fell 
lines are weak and hence the derived metallicity is practi- 
cally independent of the microturbulence. For such stars 
we have assumed £ m i cr o = 1-5 kms -1 . For the more metal- 
rich stars £ m icro has been derived by requesting that the 
derived [Fe/H] values should be independent of equiva- 
lent width. The two medium strong lines, 5197.58 A and 
5234.62 A, get high weight in this determination. The er- 
ror of £ m icro is about ±0.2 kms -1 . Due to the dominance 
of weaker lines the corresponding error of the mean value 



of [Fe/H] is negligible. It should be noted, by the way, that 
the derived microturbulence increases in a systematic way 
along the main sequence of stars in the HR-diagram, i.e. 
from about 0.8 kms -1 at T off ~ 5500 K to ~ 1.5 kms -1 
for turnoff stars at T e s ~ 6300 K. 

The error of [Fc/H] induced by the uncertainty in T e g 
is small, i.e. less than 0.02 dex. An error in logg (±0.15) 
corresponds, however, to an error of ±0.06 dex in [Fe/H]. 
Added (quadratically) to the error indicated by the line- 
to-line scatter from equivalent width measurements and 
the g /-values, we get a typical error of ±0.08 dex for 
[Fe/H] . 



4.2. Oxygen abundance from the [O i] 6300.3 A line 

The equivalent width of the 6300.3 A spectral feature was 
measured relative to the continuum defined by the cu- 
bic spline fitting procedure described in Sect. 2 and by 
integrating the line depth over an interval of 0.35A cen- 
tered on 6300.31 A. From the photon statistics we esti- 
mate an error u(W) ~ ±0.3 mA for the April spectra 
and cr(W) ~ ±0.2 mA for the July spectra. As described 
later this is confirmed by a % 2 fitting between observed 
and synthetic spectra. The total equivalent width of the 
6300.3 A feature is given in Table 4 for all stars in which 
the line is detectable (W > 0.4mA). 



1 Table Al with equivalent widths of all lines measured is 
only available in electronic form at the CDS via anonymous ftp 



to cdsarc.u-strasbg.fr (1 30.79.128.5) or via http : / /cdsweb ,u- 
3trasbg.fr / Abstract .html 



As recently shown by Allende Prieto et al. (2001), the 
[Oi] line at 6300.304 A is significantly blended by a Nil 
line at 6300.339 A with x = 4.27 eV. From a 3D model 
atmosphere synthesis of the 6 300.3 A feature in the solar 
flux spectrum of Kurucz et al. ( 1984 ) they estimate the os- 
cillator strength of the Ni line to be loggf = —2.31 for an 
adopted solar Ni abundance of loge(Ni) = 6.25. The cor- 
responding (3D) oxygen abundance of the Sun is loge(O) 
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Table 2. The derived values of effective temperature, sur- 
face gravity, metal abundance, microturbulence, absolute 
magnitude and mass. 



Dial 


Jeff 


[cgs] 


[Fe/H] 


£micro 

[km/s] 


M v 


Mass 
[M Q ] 


HD 3567 


6000 


4.07 


_2 




1 


r, 
u 


4 




n 
u 


oz 


HD 19445 


5910 


4.40 


_2 




1 


■t 
o 


r, 
O 


1 3 

J-O 


n 
u 


70 


HD 76932 


5857 


4.15 


n 


ou 




o 
z 


4 
-t 


zy 


u 


OO 


HD 97320 


5976 


4.16 


_l 


21 


1 


Q 
O 


4 


'-11 

OX 



u 


78 
I o 


HD 97916 


6328 


4.11 


_Q 


86 


I 


g 


3 


61 


1 


02 


HD 103723 


6034 


4.26 


n 
u 


70 


i 
± 


o 


4 
-t 


ou 


n 
u 


oo 


HD 106038 


5918 


4.31 


_2 


O 1 


1 


2 


4 


0'-± 


n 
u 


70 
i u 


HD 111980 


5694 


3.99 




07 


I 


2 


4 


|| 


o 


79 


HD 113679 


5678 


4.04 


u 


uu 


± 


o 
z 


1 


on 

uu 


n 

U 


00 
yu 


HD 121004 


5635 


4.34 


n 
u 


7 r i 
/ o 


± 


1 


r, 
o 


UZ 


n 

u 


SO 

OU 


HD 122196 


6107 


3.93 


_| 


70 


1 


g 


3 


66 


o 


78 


HD 126681 


5524 


4.48 


_2 


1 7 


n 
u 


7 


r, 
O 


UA 


n 
u 


70 


HD 132475 


5818 


3.95 


_2 


4:0 


| 


4 


o 


07 


n 
u 


7 r i 
1 o 


HD 140283 


5690 


3.69 


_2 


42 


| 


5 


3 


47 


o 


75 


HD 160617 


5931 


3.77 


_l 


7Q 


| 


r. 


•) 

o 


oo 



u 


89 


HD 166913 


6039 


4.11 


_l 


ou 


| 


Q 
O 


4 


22 



u 


7'-! 
* o 


HD 175179 


5643 


4.20 


_Q 


70 


| 


o 


4 


66 


o 


80 


HD 188510 


5452 


4.53 


_2 


fil 


n 
u 


a 
o 


r, 
u 


0'-± 


n 
u 


R8 


HD 189558 


5613 


3.91 


_2 


|| 


| 


2 


4 


03 

Uu 


n 

u 


70 


HD 195633 


5987 


3.95 


_Q 


54 


| 


4 


3 


38 


1 


10 


HD 205650 


5733 


4.39 


_l 


1 7 


| 


n 
u 


r, 
o 


21 



u 


70 


HD 213657 


6114 


3.85 


— 1 


Q 1 


1 
1 


c 



Q 

o 


1 o 


U 




HD 298986 


6071 


4.21 


-1 


35 


1 


3 


4 


41 





76 


HD 338529 


6239 


3.87 


-2 


34 


1 


5 


3 


40 





79 


CD -30° 18140 


6274 


4.12 


-1 


91 


1 


5 


4 


06 





75 


CD -35° 14849 


6117 


4.12 


-2 


34 


1 


5 


4 


24 





70 


CD —57° 1633 


5909 


4.32 


-0 


89 


1 


2 


4 


66 





84 


G 013-009 


6338 


4.00 


-2 


30 


1 


5 


3 


70 





76 


G 020-024 


6440 


4.30 


-1 


78 


1 


5 


4 


33 





78 


G 075-031 


5952 


4.07 


-1 


04 


1 


3 


3 


96 





87 


G 088-032 


6251 


3.81 


-2 


41 


1 


5 


3 


26 





80 


G 126-062 


5941 


3.98 


-1 


58 


1 


3 


3 


93 





76 


G 183-011 


6092 


4.04 


-2 


10 


1 


5 


4 


08 





70 


G 271-162 


6230 


3.93 


-2 


31 


1 


5 


3 


61 





75 


LP 815-43 


6514 


4.23 


-2 


72 


1 


5 


4 


17 





75 



= 8.69 when using the accurate oscillator strength of the 
[O I] line, loggf = 



-9.72 (Allcndc Prieto et al. |200l| ). 
The 6300.3 A feature in the solar flux spectrum has 
an equivalent width of 5.5 mA. Subtracting the calculated 
equivalent width of the Ni line (1.4 mA), we determine 
a (ID) oxygen abundance of the Sun, loge(O) = 8.74, in 
quite good agreement with the 3D result of Allcnde Prieto 
et al., but substantially lower than previous determina- 
tions that did not take the Ni line into account. Indeed, 
when accounting for the difference of — 0.04dex between 
3D and MARCS model atmosphere for the solar O abun- 
dance, our analysis agrees almost perfectly with that of 
Allende Prieto et al. ( [200 

Adopting now \oggf = —2.31 for the Nil line and as- 
suming that the Ni abundance scal es pr oportional to the 
Fe abundance (Edvardsson et al. 1993 ), the equivalent 



Table 4. Equivalent widths, W to t, of the 6300.3 A absorp- 
tion feature together with calculated values for the blend- 
ing Nil line, Wm, and values of loge(O) and [O/Fe] as de- 
rived from ID model atmospheres. The last two columns 
give the ±lcr errors of [O/Fc]. 



Star [Fe/H] 



Wtot 

mA 



mA 



loge(O) [O/Fe] +a 



Sun (flux) 





00 


5 


5 


1 


4 


8 


74 





00 








TTT\ or/ 1 

HD 3567 


— 1 


16 


1 


7 





1 


8 


00 





42 


+0 


08 


—0.10 


HD 76932 


-0 


86 


3 


3 





2 


8 


32 





44 


+0 


06 


—0.06 


1 1 1 J y / ozu 


— 1 


9 1 


i 
i 


9 
Z 





i 
i 


7 
/ 


oo 


u 


Q 

OO 


n 


1 1 
1 1 


n ic 
— U.IO 


HD 103723 


-0 


79 


1 


7 





2 


8 


13 





18 


+0 


09 


-0.10 


HD 111980 


-1 


07 


2 


1 





1 


7 


96 





29 


+0 


08 


-0.09 


HD 113679 


-0 


66 


4 


5 





4 


8 


39 





31 


+0 


05 


-0.05 


HD 121004 


-0 


75 


3 


1 





3 


8 


30 





31 


+0 


06 


-0.07 


HD 126681 


-1 


17 


1 


8 





1 


8 


03 





46 


+0 


09 


-0.10 


HD 132475 


-1 


45 


1 


5 





1 


7 


78 





49 


+0 


09 


-0.11 


HD 140283 


-2 


42 





5 








7 


09 





77 


+0 


15 


-0.22 


HD 160617 


-1 


79 





6 








7 


37 





42 


+0 


19 


-0.30 


HD 166913 


-1 


56 





6 








7 


57 





39 


+0 


18 


-0.30 


HD 189558 


-1 


11 


3 








1 


7 


06 





43 


+0 


06 


-0.07 


HD 195633 


-0 


54 


3 








4 


8 


28 





08 


+0 


07 


-0.07 


HD 205650 


-1 


17 


1 


6 





1 


8 


00 





43 


+0 


09 


-0.11 


HD 213657 


-1 


94 





4 








7 


31 





51 


+0 


18 


-0.30 


HD 298986 


-1 


35 





8 








7 


77 





38 


+0 


15 


-0.21 


G 075-031 


-1 


04 


1 


2 





1 


7 


85 





15 


+0 


11 


-0.15 



width of the Ni line (see Table 4) has been calculated and 



subtracted from the total equivalent width of the 6300.3 
A feature. The resulting VF([0 1]) is used to determine the 
oxygen abundances of the stars relative to that of the Sun, 
i.e. [O/H]. As seen from the table, the Ni line is rather 
weak in the metal poor stars and its inclusion decreases 
the derived oxygen abundance by less than 0.04 dex for 
stars with [Fe/H] < —1.0. For the Sun, however, the effect 
of including the Ni I line is a decrease of the oxygen abun- 
dance by 0.13 dex. Hence, the net effect of the Ni line for 
metal-poor halo stars is an increase of [O/Fe] with about 
0.1 dex. 

The effect of uncertainties in T e g and logg on [O/Fe] 
is small; cr(T G ff) = ±70 K corresponds to er([0/Fe]) = 
±0.035 dex, and most important, the effect of varying the 
gravity is very nearly the same on loge(O) and loge(Fe), 
when the two abundances are derived from the [O i] line 
and Fe n lines, respectively. For example, an increase of 
the gravity of HD 132475 by as much as 0.3 dex changes 
[O/Fe] by less than 0.01 dex. Hence, [O/Fe] is not affected 
by errors in logg. In practice, the error budget of [O/Fe] is 
totally dominated by the error of the equivalent width of 
the [O i] line. In order to check our estimate of this error 
source we have made a \ 2 analysis of a comparison of syn- 
thetic spectra with some representative observed spectra. 

The synthetic spectra are computed for the spectral re- 
gion 6299 - 6302 A. In addition to the [Oi] and Nil lines, 
Si I 6299.58 A, Sen 6300.68 A, and Fei 6301.50 A are in- 
cluded, see Fig. |i[ The loggf values of these lines were set 
to get agreement between their calculated and observed 
equivalent widths, but they have practically no influence 
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Fig. 4. ID model atmosphere synthesis of the observed 
spectrum of HD 3567 {S/N ~ 540). The xL« is obtained 
for loge(O) = 7.95 (full drawn line). The two dashed lines 
correspond to oxygen abundances changed with ±0.2 dex. 
The line and continuum regions used for the x 2 analysis 
are indicated. Note, that the two deviating points to the 
right of the Fe I line are due to telluric [O i] emission. 



on the x 2 analysis, because the spectral regions of these 
lines are not included. The synthetic spectra are folded 
with a Gaussian profile representing the instrumental, 
macroturbulent and rotational broadening. The FWHM 
of this Gaussian is determined from the Fel 6301.50 A 
line, e.g. 5.8 kms -1 in the case of HD 3567. The x 2 analy- 
sis is then performed over the spectral regions indicated on 
Fig. 4, i.e. a total of 42 datapoints, of which the C regions 
control the continuum setting and the L region determines 
the O abundance. For each chosen oxygen abundance the 
continuum setting is optimized to get the smallest possible 
X 2 ■ In the case of HD 3567 the minimum x 2 1S obtained 
for loge(O) = 7.95 ±0.05. The reduced x 2 is 0.88, and the 
quoted (one-sigma) error corresponds to Ax 2 = 1. This 
oxygen abundance is 0.05 dex lower than the value deter- 
mined from the measured equivalent width of the 6300.3 
A feature, but well within the estimated error bars. A sim- 
ilar x 2 analysis was performed for HD 160617 (July 2000 
spectrum) and HD 132475 and HD 189558 (April 2000) 
confirming in all cases the oxygen abundances found from 
the measured equivalent widths of the 6300.3 A spectral 
feature. 

4.3. Oxygen abundance from the 1 triplet 

For the 15 stars belonging to the July 2000 program the O I 
triplet at 7774 A was included in the UVES spectra. The 
triplet occurs in two consecutive echelle orders each with 
a S/N of 200-300 in the wavelength region of the triplet. 
Hence, two sets of equivalent widths could be measured 
and from the comparison, the error of the mean equiva- 
lent widths (given in Table 5) is estimated to be less than 
1 mA. The high quality of the spectra is evident from Fig. 




7776 



7770 



7772 7774 
Wavelength (A) 



7776 



Fig. 5. Spectra of two stars in the Oi triplet region 



S, which shows the average spectra of the triplet for two 
representative stars. 

From the equivalent widths of the triplet, LTE abun- 
dances were first derived using the model atmospheres 
employed in the [Oi] analysis. The adopt ed ^/-values are 
those recommended by Wiese et al. ( 1996 ): log g/ = 0.369, 
0.223, and 0.002. For a given star the three Oi lines pro- 
vide the same oxygen abundance to within about 0.03 
dex except in the case of LP 815-43, for which the disper- 
sion is around 0.09 dex due to the faintness of the lines. 
Estimated uncertainties of ± 70 K in T c g, and ± 0.15 in 
log 5 correspond to Aloge(O) of ± 0.05 and ± 0.06, re- 
spectively. Uncertainties arising from the microturbulcncc 
and [Fe/H] may be neglected. We estimate the total er- 
ror of log e(O) to be about ± 0.1 dex, whereas the error of 
[O/Fe] is in fact a bit smaller, because the gravity induced 
errors on the O and Fe abundances nearly cancel. 



It is well established that the line formation process 
of the O I triplet in late-type stars is affected by depar- 
tures from LTE (e.g. Kiselman 2001). To account for these 
important non-LTE effects, we have performed statisti- 
cal equilibrium calculations using version 2.2 of the pub- 
licly available code MULTI (Carlsson 1986) for the stars 
listed in Table 5. For the purpose, a 23 level model atom 
(22 levels for O I plus the ground state of O II) with 43 
bound-bound and 22 bound-free radiative transitions has 
been adopted based on a model atom kindly provided by 
Dan Kiselman (2000, private communication). The radia- 
tive data comes mainly from Opacity Project calculations 
(Seaton 1987). As evidence is mo unting that the simple- 
minded recipe of Drawin ( 1968 ) to incorporate the ef- 
fects of inelastic collisions with hydrogen overestimates 
the cross- sectio ns by at least two orders of magn itude 
(Lambert |l993| ; Belyayev et al. |l999j ; Kiselman |200l| ; Paul 
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Table 5. Equivalent widths of the 7771.2, 7774.2, and 
7775.4 A lines in the O I triplet and the derived mean LTE 
and non-LTE oxygen abundances. The last column gives 
[O/Fe] based on the non-LTE oxygen abundance. 



Star [Fe/H] 



Wi 
mA 



U. 

mA 



W 3 loge(O) loge(O) 
mA 



LTE 



n-LTE 



[O/Fe] 

n-LTE 



-L1J-/ OOvJ 1 
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1 f\ 
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41 
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90 
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14 
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n 
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HD 19445 


-2 


04 


13 


9 


9 


7 


6 


3 


7 


49 


7 


39 





69 


HD 106038 


-1 


37 


39 


5 


32 


1 


22 


9 


8 


10 


7 


94 





57 


HD 140283 


-2 


42 


7 


1 


5 


8 


3 


2 


7 


11 


7 


02 





70 


HD 160617 


-1 


79 


19 


9 


14 


9 


9 


7 


7 


46 


7 


34 





39 


HD 213657 


-1 


94 


26 


5 


20 


7 


14 





7 


56 


7 


41 





61 


HD 338529 


-2 


34 


16 


7 


13 


3 


7 


4 


7 


20 


7 


06 





66 


-30° 18140 


-1 


91 


27 


2 


21 


5 


14 


7 


7 


56 


7 


41 





58 


-35° 14849 


-2 


34 


11 


1 


8 





5 


2 


7 


15 


7 


04 





64 


G 013-009 


-2 


30 


14 


8 


12 


4 


7 


8 


7 


15 


7 


01 





57 


G 020-024 


-1 


78 


26 


8 


20 


6 


13 


6 


7 


49 


7 


34 





38 


G 075-031 


-1 


04 


55 


7 


44 


4 


34 


3 


8 


28 


8 


06 





36 


G 126-062 


-1 


58 


38 


4 


30 


1 


21 





7 


96 


7 


79 





63 


G 271-162 


-2 


31 


13 


2 


9 


5 


5 


9 


7 


08 


6 


96 





53 


LP 815-43 


-2 


72 


6 





4 


7 


3 


9 


6 


71 


6 


46 





43 



Barklem, 2001, private communication), such collisions 
are not included in the calculations. We have, however, 
performed test calculations with H collisions estimated fol- 
lowing Drawin's recipe multiplied with different factors for 
all stars, as discussed briefly below. Collisional excitation 
and ionization with electrons and charge transfer collisions 
with hydrogen and protons from the ground levels of O I 
and O II are included. For details of the adopted atomic 
data w e refe r to Kiselman (1991, 1993) and Carlsson & 
Judge ( 1993 ). We note that the non-LTE results for O I are 
quite insensitive to including additional levels and transi- 
tions as the O I triplet line formation is well described by 
a two-level atom line source function with the line opacity 
remaining close to its LTE value (Kiselman |l993| ). The 
important factors for the line formation instead is the line 
itself: the mean intensity dropping below the Planck func- 
tion in the infrared, the collisional cross-sections and the 
line strength. 

The ID LTE abundances for the three O I triplet lines 
have been individually corrected for the predicted depar- 
tures from LTE and the resulting mean non-LTE oxygen 
abundance and the corresponding [O/Fe] value are given 
in the last two columns of Table 5. In all cases the non- 
LTE effects have been estimated using MARCS model at- 



mospheres (Asplund et al. 1997) with the stellar parame- 
ters listed in Table 2. The non-LTE abundance corrections 
amount to a 0.1-0.25 dex decrease of the LTE results, in 
agreement with previous calculations without large H col- 
lisional cross-sections (e.g. Kiselman 1991, 1993). The re- 
sults remain largely intact also when including the effects 
of H collisions. With a factor of 0.01 times the Drawin esti- 
mates, the non-LTE abundance corrections are within 0.01 
dex of the case without H collisions. Even when applying 
the Drawin values directly, the corrections only change by 
< 0.04 dex with the exception of LP 815-43 where the dif- 



ference amount to 0.1 dex. We note that LP 815-43 has the 
largest non-LTE effect (—0.25 dex), which becomes even 
larger at yet smaller metallicities. 

4.4. The influence of stellar granulation on [Oi] and 
Fen lines 

Although the [O i] and Fe II lines are virtually immune to 
departures from LTE, the strengths of the lines depend on 
the detailed structure of the adopted atmospheric model. 
Properties of particular importance are the temperature 
gradient, the ratio of the line and continuous opacity and 
the atmospheric inhomogcncitics and velocity fields intro- 
duced by convection in the surface layers, also known as 
stellar granulation. Since the [O i] and Fe II lines are sen- 
sitive to different parts of the atmospheres due to the dif- 
ference in excitation potential and ionization balance, it 
can not be taken for granted that the model atmosphere 
dependence cancels out when forming [O/Fe] abundance 
ratios, as is the case for analyses using ID model atmo- 
spheres. 

In order to estimate the effects of stellar granulation on 
our derived O and Fe abundances, the new generation of 
3D, time-dependent hydrodynamical model atmospheres 
(e.g. Asplund et al. 1999, 2000a,b,c) has been applied to 
study the line formation of the [Oi] and Fell lines. The 
procedure is the same as that adopted for our recent in- 
vestigation of the OH line formation (Asplund & Garcia 
Perez 2001): a strictly differential comparison has been 



performed between ID MARCS and 3D model atmospheres 
with identical stellar parameters. For details of the con- 
vection simula tions the reader is referred to Asplund & 
Garcia Perez (2001). The models employed here have a 
range of metallicities ([Fe/H] = 0.0 to —2.0 and corre- 
spond to solar-like stars (T c s ~ 5800, logg = 4.4) and 
turn-off stars (T c ff ~ 6200, log g = 4.0). In addition 
we have included the subgiant HD 140283 (T cff ~ 5700, 
logg = 3.7, [Fe/H] = —2.5 ) usin g the 3D simulation 
presented in Asplund et al. ( 1999 ). The same spectrum 
synthesis and atomic data (transition properties, continu- 
ous opacities, equation-of-state and chemical composition) 
have been employed in both ID and 3D to isolate the im- 
pact of convection treatment from other possible sources 
of uncertainties. For all line calculations the assumption 
of LTE has been made, which, as noted before, is well jus- 
tified. In the 3D case no macro- or microturbulence are 
needed while £ m i cr o = 1.0 km s" 1 has been adopted for the 
corresponding ID calculations although only of some im- 
portance for the Fen lines in the solar metallicity models. 

The results for th e [O i] line and a subset of the Fell 
lines utilized in Sect. 4T are summarized in Table ^. The 
impact on [O /Fe] is presented in Fig. ^| where the mean 
3D granulation corrections for the different T e g models are 
shown relative to the solar case 0. Even if the influence of 
granulation on the derived O and Fe abundances are rela- 



2 The r esults shown in Fig. ^| supersede those presented in 
Asplund (|200lj ) by the inclusion of Fell lines, use of addi- 
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Table 6. Comparison of the oxygen and iron LTE abun- 
dances derived with ID hydrostatic and 3D hydrodynam- 
ical model atmospheres. The 3D abundances are those 
which reproduce the equivalent widths of the [O i] and 
Fe II lines computed using a ID model atmosphere, a mi- 
croturbulence of £ m i cro = 1.0 km s -1 and the abundances 
given in the third and fourth columns. 







[Fe/H] 


[O/Fef 


Alog £ (0) 


Aloge(Fe) b 


[K] 


[cgs] 










5767 


4.44 


+0.0 


±0.00 
+0.20 
+0.40 


-0.04 
-0.04 
-0.03 


-0.02 


5822 


4.44 


-1.0 


+0.20 
+0.40 
+0.60 


-0.08 
-0.08 
-0.06 


+0.04 


5837 


4.44 


-2.0 


+0.40 
+0.60 
+0.80 


-0.12 
-0.12 
-0.12 


+0.10 


6191 


4.04 


+0.0 


+0.00 
+0.20 
+0.40 


-0.01 
-0.01 
-0.01 


-0.01 


6180 


4.04 


-1.0 


+0.20 
+0.40 
+0.60 


-0.07 
-0.11 
-0.11 


+0.02 


6178 


4.04 


-2.0 


+0.40 
+0.60 
+0.80 


-0.16 
-0.17 
-0.15 


+0.06 


5690 


3.67 


-2.5 


+0.40 
+0.60 
+0.80 


-0.22 
-0.24 
-0.24 


+0.05 



For these calculations loge(O) = 8.70 and loge(Fe) = 7.50 
have been adopted for the Sun. 

The impact of 3D model atmospheres on the derived stel- 
lar metallicities have been investigated for in total six Fe II 
lines (6149, 6238, 6247, 6417, 6432 and 6456 A). In all cases 
the different Fell lines give the same granulation correc- 
tions to within 0.03 dex. 



tively small in most cases (~ 0.1 dex), the combined effects 
on [O/Fe] can be substantial amounting to ~ —0.2 dex at 
[Fe/H] = —2.0 due to the opposite behaviour of the [Oi] 
and Fell lines. In particular, the 3D correction to [O/Fc] 
is ~ -0.26 dex for HD 140283. 

The selected Fen lines are predominantly sensitive to 
the temperature gradient in the deep atmospheric lay- 
ers making them somewhat weaker in 3D at low [Fe/H] 
(Asplund et al. 



1999) 



The [Oi] line on the other hand 
feels the cool, sub-radiative equilibrium upper layers at 
low metallicities, making it significantly stronger in 3D 
than in ID model atmospheres. At first sight, this be- 
haviour may seem peculiar since the population of the 
lower level of the [O i] line is essentially insensitive to the 
temperature for these stars because it originates from the 
ground state of the dominant ionization stage. Instead the 
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tional and more temporally extended 3D model atmospheres 
and shown relative to the solar case. 



Fig. 6. The effects of stellar granulation on [O/Fe] based 
on [O i] and Fe II lines according to 3D hydrodynamical 
model atmospheres. Due to the opposite behaviour of the 
[O i] and Fe II lines the impact on [O /Fe] can be substantial 
at low metallicities. For these calculations the adopted so- 
lar O and Fe abundances were 8.70 and 7.50, respectively. 

main effect on the line comes not from a change in line 
opacity but from differences in continuous opacities be- 
tween 3D and ID models. The strengths of weak lines like 
the [Oi] (and Fen) lines depend on W\ oc 4 me / K ™ nt 00 
Noilx = eV)/N m N e ~ N /N H -1/N C since the continu- 
ous opacity is provided mainly by H~ around 630 nm and 
oxygen and hydrogen are predominantly in neutral form 
for these atmospheres. At low metallicities and the rele- 
vant temperatures, the number density of free electrons 
is very temperature sensitive. Thus, in order to compen- 
sate for the lack of electrons in the cool surface layers in 
the 3D model atmospheres, the abundance No/N-a must 
be lowered to retain the same equivalent width as in the 
corresponding ID model (Fig. Q). 

This effect is not restricted to the [O i] line but applies 
equally well to any resonance line of a majority species, 
as illustrated in Fig. |^. For this comparison a number of 
fictitious lines of Oi, Fell and Si with different excita- 
tion potentials have been computed for a few snapshots in 
the HD 140283 simulation as well as for the correspond- 
ing ID MARCS model atmosphere. The ^/-values for the 
different transitions have been adjusted to produce a line 
of similar strength as the real [Oi] line, i.e. 0.5mA. In all 
cases, the lines of Oi, Fen and Si behave almost iden- 
tically since all are majority species. For low excitation 
potentials the effect on the continuous opacities dominate 
while higher excitation lines, like the Fe n lines utilized in 
the present study, are more sensitive to the effects on the 
line opacities and thus the temperature gradient in the 
deeper atmospheric layers. 

For the discussion below in terms of the true mean 
trend in [O/Fe] as a function of [Fe/H] it is reasonable to 
approximate the 3D results according to the simple rela- 
tion A [O/Fe] ~ 0.11 • [Fe/H]. We emphasize that this rela- 
tion can not be assumed to be automatically valid for stel- 
lar parameters differing greatly from those given in Table 
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Fig. 7. Temperature (upper panel) and electron pressure 
(lower panel) as a function of continuum optical depth at 
500 nm in one snapshot of the 3D hydrodynamical sim- 
ulation of HD 140283. Also shown is the corresponding 
ID MARCS photospheric structure (crosses connected with 
solid lines), illustrating the much lower temperatures and 
consequently electron pressures encountered in 3D model 
atmospheres in the upper atmospheric layers. 
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Fig. 8. The 3D LTE abundance corrections (3D- ID) for 
a selection of fictitious O I, Fe II and S I lines, all with an 
equivalent width of 0.5 mA and located at 630 nm but with 
different excitation potentials, for the metal-poor subgiant 
HD 140283 ([Fe/H]= -2.5). The lines of these majority 
species behave almost identically with low excitation po- 
tential lines being significantly stronger in 3D than in ID 
model atmospheres. 



q. In particular, giants and subgiants, like BD +23°3130 
which is discussed at some length below, may diverge from 
this relation given their lower T e g and log g. The sign is 
unlikely to change, however, and therefore [O/Fe] ratios 
based on [O i] and Fe II lines should always become smaller 
in 3D than in ID. 



5. Discussion 

5.1. [Oi] and the oxygen problem 

Our derived [O/Fe] values based upon the ID hydrostatic 
model atmospheres are plotted versus [Fe/H] in Fi g. p| to - 
gether with disk stars from Nissen & Edvardsson ( 1992 ). 
Their [O /Fe] values were also obtained from the [O i] 
6300.3 A and Fell lines, but without any correction for 
the effect of the blending Nil line. Hence, we calculated 
the equivalent width of the blend, subtracted it from the 
observed equivalent width of the 6300.3 A feature, and re- 
calculated the oxygen abundance. It is these new [O/Fc] 
values for the disk stars that are plotted in Fig. |[ No 
attempt has been made to correct for the small offset be- 
tween the T e ff scale of Alonso et al. used here and that used 
by Nissen & Edvardsson, which was based on Edvardsson 
et al.'s (1993) theoretical calibration of the b 



y index. 

Although an offset between the two sets of [O/Fe] values 
may be present, we expect it to be small due to the low 
sensitivity of [O/Fe] to T e g. 

We emphasise again that the combination of [O i] and 
Fen lines, as analysed with ID model atmospheres yields a 
result for [O/Fe] which is insensitive to reasonable errors 
in the atmospheric parameters and to non-LTE effects. 
The [Fe/H] derived from the Fell lines is sensitive to the 
adopted surface gravity but insensitive to non-LTE effects. 
Our surface gravity estimates are not determined spec- 
troscopically, as are many published values using the as- 
sumption of ionisation equilibrium. Although other spec- 
troscopic indicators of the oxygen (and iron) abundances 
may be affected by multiple sources of uncertainty, the 
leading source of a systematic error in our abundance de- 
terminations is likely to arise from the simplifying assump- 
tions entering into the calculation of the ID model atmo- 
spheres. 

As discussed in Sect. 4.4 and highlighted in Fig. [| 
the effect of the stellar granulation as represented by the 
3D hydrodynamic model atmospheres is significant for the 
most metal-poor stars. While emphasising that the grid 
of 3D models is as yet incomplete, we have interpolated 
corrections from Table |^ and constructed Fig. [l^ as an ap- 
proximation to the [O/Fe] versus [Fe/H] for the 3D model 
atmospheres. 

Our sample is the first using accurate measurements of 
the [O i] line to provide extensive results for a large number 
of metal-poor (say, [Fe/H] < — 1) dwarfs and subgiants. In 
light of the fact that published results using the [O i] line 
are mostly based on inferior spectra, a detailed comp arison 
is not attempted; for example, Cavallo et al. ( |l997| ) draw 
from the literature two quite different measurements for 
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the 6300 A line in HD 76932 (W = 10 mA, and 1.8 mA) 
with neither one close to our measurement of 3.3 mA. Two 
metal-poor subgiants were analysed recently by Fulbright 
& Kraft ( |1999|) . Since one - BD+23°3130 - has been the 
object of several recent discussions, we discuss it in some 
detail later in this section. Results from the [Oi] lines in 
giants are discussed separately. 

Attention was sharply focussed on oxygen abundances 
in metal-poor stars by studies of the OH ultraviolet lines 
(Israelian et al. |199S| , |200l| ; Boesgaard et al. |1999| ). These 
showed that [O/Fe] indices derived from the OH ultra- 
violet lines in metal-poor dwarfs were higher than the 
values obtained previously from [O i] lines in metal-poor 
giants. The latter results had shown that [O/Fe] was ap- 
proximately constant (~ 0.4) for [Fe/H] < —1. In sharp 
contrast, the OH-based results showed [O/Fe] to increase 
quasi- linearly with decreasing [Fe/H]. The discrepancy be- 
tween the abundance from OH in dwarfs and [O i] in giants 
was negligible for [Fe/H] > — 1 but increased steadily to 
about 0.6 to 0.7 dex at [Fe/H] ~ —3. Earlier, an analysis 
of the Oi triplet lines near 7774 A in metal-poor dwarfs 
had also given high [O/Fe] values indicating a sharp dis- 
agreement with result s from the [O i] 6300 A line in giants 
(Abia fc Rebolo |l989| ). 

With our analyses of the [O i] line in dwarfs, we look 
afresh at the seemingly discrepant oxygen abundances. In 
particular, we focus on three key questions: 



— Do the [O i] line and the OH lines provide consistent 
oxygen abundances for dwarfs and subgiants? 

— Do the [Oi] line and the Oi 7774 A triplet provide 
similar abundances for dwarfs and subgiants? 

— Is the oxygen abundance - metallicity dependence sim- 
ilar for dwarfs and giants? 

5.2. [O i] and OH lines 

How do our new [O i] -based abundances agree with those 
derived from the OH ultraviolet lines, when both are ob- 



Fig. 9. [O/Fe] vs. [Fe/H] as derived on 
the basis of ID model atmospheres. 
Filled circles refer to stars from the 
present paper and open circles t o dis k 
stars from Nissen & Edvardsson (1992), 
in both cases with the oxygen abun- 
dance determined from the [O i] 6300.3 
A line. The triangles refer to halo stars 
with non-LTE oxygen abundances de- 
termined from the O I triplet at 7774 A. 



tained from classical ID model atmospheres? A compar- 
ison is illuminating given the heat the 'oxygen problem' 
has generated in recent times. 

Eight of our programme stars were analysed by 
Israelian et al. ( 1998j ) or Israelian et al. (2001). Differences 
in adopted model atmosphere parameters are small: 
AT off = 91± 103 K, A log g = 0.18 ± 0.26, and A [Fe/H] 
= 0.19 ± 0.07 where the sense of the difference A is 'this 
paper — Israelian et al.'. Israelian et al. chose their model 
atmosphere parameters and [Fe/H] from the literature. 
Different grids of ID model atmospheres were used in the 
two studies. 



For four of the eight common stars, we determined 
the oxygen abundance from the [O i] line. For this quartet 
which includes HD 140283 - the most metal-poor star in 
our sample - the average difference between our abun- 
dances ([O/H] not [O/Fe]) from the ID model atmo- 
spheres (Table 4) and those from OH is —0.06 dex. If ad- 
justments are made for the use of different grids of model 
atmospheres and atmospheric parameters, the mean dif- 
ference may increase to about —0.10 dex for these stars 
with [O/H] between -0.4 and -1.6. Israelian et al. Ql998|) 



adjusted the ^/-values of the OH lines to match the solar 
OH lines. In analysing spectra of the metal-poor stars, it 
seems preferable to adopt not solar g /-values but the re- 
l iablc theoretical-experimental g /-values. Melendez et al. 
( 2001 ) argue that were the latter gf- values adopted the 
published OH-based abundances would be lowered by 0.1 
- 0.2 dex, and then the mean difference would increase 
from —0.1 dex to between 0.0 and 0.1 dex. 

Six of our stars are in the sample analysed by 
Boesgaard et al. ( 1999] ), who presented OH-based results 
for two different T e g scales. Their atmospheric parame- 
ters based on Carney's ( |1983b[ ) scale are similar to ours: 
AT off = 15 ± 34 K, A log 5 = 0.09 ± 0.20, and A [Fe/H] 
= 0.12 ± 0.07. Of these six, we determined oxygen abun- 
dances for two (including HD 140283). After small cor- 
rections are applied for differences in the adopted model 
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atmospheres, differences in [O/H] are +0.1 and —0.1 dex. 
Since Boesgaard et al. also adjusted ^/-values to fit so- 
lar lines, there will be, as in the case of Israelian et al., 
a slight revision downward of the published abundances 
when laboratory (//-values are used. ^] 

These two comparisons show that, where a direct com- 
parison is possible, analyses of the ultraviolet OH and the 
6300 A [O i] lines using 1-D model atmospheres give con- 
sistent oxygen abundances for metal-poor dwarfs and sub- 
giants. This result applies to stars with [Fe/H] > —2.5. 

Comparison of results is more commonly made in 
terms of the run of [O/Fe] versus [Fe/H], but this intro- 
duces other considerations. The adopted solar O and Fe 
abundances now enter the picture. Our adopted solar oxy- 
gen abundance is about 0.2 dex smaller than chosen in the 
papers covering OH. Since the OH analyses used solar gf- 
values, the lower solar abundance should not affect their 
[O/Fe] values but it likely accounts for the fact that the 
solar gj-values for OH lines were 0.1 - 0.2 dex smaller 
than the laboratory (//-values (see above). Additional ef- 
fects should be noted. Kurucz models were used for the 
OH analyses but MARCS models for the [Oi] line, but 
the effect on [O/Fe] indices is likely small (< 0.1 dex). 
Also, our [Fe/H] are systematically slightly larger than 
those adopted in the OH analyses. The net effect of these 
various factors may not exceed about 0.1 dex in [O/Fc]. 

The OH-based results imply a linear relation between 
[O /Fe] a nd [Fe/H] with a slope of -0.33± 0.02 (Israelian et 
al. pOOl ), a relation which runs along the upper boundary 
of the points in Fig. ^for [Fe/H] < —1.5. In short, the two 
indicators give reasonably consistent results. 

The much discussed subgiant star BD +23°3130, which 
has been added to Fig. 0, was thrust into the oxygen de- 



bate by Fulbright & Kraft (1999) who measured the [O i] 
6300 A line and argued that the resulting [O/Fe] did not 
support the high values from OH ultraviolet lines. Others 
have discussed the [Oi] line in BD+23°3130 obtaining 
different results. We have reanalyzed BD +23°3130 in the 
same way as our 'own' stars. We take the eq uivalent width 
of the 6300 A line from Cayrel et al. ( 2001 ) who obtained 
an UVES spectrum of BD +23°3130 and determined the 
equivalent width of [Oi] 6300.3 A to W = 1.5 ± 0.5mA, 
where the quoted ler error seems to be a rather conserva- 
tive estimate in view of the published high S/N spectrum. 

From unpublished Stromgren photometry of 
BD+23°3130 (Schuster et al. pOO^ ) T cS (b - y) = 5195 K 
is derived, and V - K = 1.975 (Laird et al. |l98SD leads 
to T cff (V - A") = 5145K. In both cases the IRFM T cB 
calibration for giants by Alonso et al. (|l99S|) has been 



3 The O abundance derived from OH UV and IR lines is 
dependent on the assumed dissociation energy of OH. Most 
analyses appear to ad opt D p = 4.392 eV, a recommendation 
by Huber & Herzber g (|l979| ) based on a spectroscopic estimate 
by Carlone & Dalby ( |1969| ). A recent r ecomm endation provides 
Do = 4.413 ± 0.003 eV(Ruscic et al. |2002j ) from new experi- 
mental and theoretical evidence and a reevaluation of the spec- 
troscopic data. Use of the new Do will result in slightly lower 
O abundances. This revision is ignored here. 



used and it was assumed that the star is unreddened. 
Adopting the average T c g — 5170 K, the gravity derived 
via the Hipparcos parallax is logg = 3.0 ± 0.25. Two 
of the Fen lines from the list in Table 3, 5197.58 A 
and 5234.63 A, were observed by Fulbright & Kraft 
( 1999Q and have equivalent widths of 19.6 and 23.0 mA, 
respectively. The corresponding metal abundance is 
[Fe/H] = —2.39. With these parameters and using the 
equivalent width of the [O i] 6300.3 A line given by Cayrel 
et al., we get [O/Fe] = 0.64 ± 0.15, the value plotted in 
Fig. |. However, it cannot be excluded that BD +23°3130 
(distance ~ 230 pc) is significantly reddened. Act ually , 
the observed Hp index of the star (Schuster et al. 2002 ) 
suggests T e ff ~ 5270 K based again on the calibration 
of Alonso et al. ( 1999 ). This temperature would lead to 
[O/Fe] = 0.69 ±0.15. 

This shows that BD +23°3130 confirms our result for 
HD 140283. The oxygen abundance for BD +23°3130 is 
loge(O) = 6.99. Israelian et al. (2001) for a Kurucz model 
having almost identical atmospheric parameters obtained 
loge(O) = 7.10 from OH ultraviolet lines. Balachandran 
et al. (2001) also adopted Cayrel et al.'s W for the 6300 
A line and with atmospheric parameters essentially iden- 
tical to ours found loge(O) = 7.03 ± 0.15. Reanalysis of 
Israelian et al.'s (1998) equivalent widths of OH ultravi- 
olet lines with laboratory (//-values gave loge(O) = 7.06 
± 0.11, and analysis of their own measurements of OH in- 
frared vibration-rotation lines gave loge(O) — 7.00 ± 0.07. 
Thus, as for other metal-poor stars, the [O i] and OH ul- 
traviolet (and infrared) lines give concordant results. 

The OH molecule also provides vibration-rotation lines 
from the ground electronic state. Fundamental or first 
overtone lines have been detected and analysed in a few 
cool metal-poor dwarfs. Given spectra of adequate S/N, 
the infrared lines have the advantages over the ultraviolet 
lines that they lie in a cleaner spectral interval, and are 
likely to be formed in LTE. In addition to BD +23°3130 
mentioned above, results are available for three additional 
metal-poor dwarfs. Balachandran & Carney ( 1996] ) discuss 
HD 103095 (Gmbl830) obtaining loge(O) = 7.83 from 
OH infrared line s or [ O/Fe] = 0.31 for [Fe/H] = -1.22. 
Melendez et al. ( E00l[ ) give a similar result: after correc- 
tion for our lower solar oxygen abundance, [O/Fe] ~ 0.4 
at [Fe/H] = —1.36. These values fit the results in Fig. || 
for [Fe/H] ~ -1.3. Results for two other dwarfs at [Fe/H] 
= — 0.8 and —1.8 according to Melendez et al. fall at the 
upper end of the distribution in Fig. ^|. Balachandran et 
al.'s run of [O/Fe] versus [Fe/H] from OH infrared lines in 
nine stars with [Fe/H] from —1 to —3 resembles the dis- 
tribution of points in Fig. but a 0.1 dex increase for the 
lower solar oxygen abundance should be applied. A closer 
comparison is not possible because details of the infrared 
analyses are unpublished. 

The preceding comparisons between [O i] and OH lines 
were based on the use of ID model atmospheres and ig- 
nored the effects of granulation. In Section 4.3, we dis- 
cussed these effects on the [Oi] line. Abundances cor- 
rected approximately for these 3D effects are shown in 
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Fig. 10. [O/Fe] vs. [Fe/H] after approx- 
imate correction for the effects of stellar 
granulation. Filled circles refer to stars 
from the present paper and open circles 
t o disk stars from Nissen & Edvardsson 
in both cases with the oxygen 
abundance determined from the [O i] 
6300.3 A line. 



Fig. [T(J. Analogous effects on the OH lines were discussed 
by Asplund & Garcia Perez (2001). For a given effective 
temperature, the infrared and ultraviolet OH lines are sim- 
ilarly sensitive to the effects of granulation (Asplund & 
Garcia Perez 2001). The effects of granulation on the OH 
lines appears to decline with decreasing effective temper- 
ature. Since stars observed for their infrared OH lines are 
generally cooler than those observed for the ultraviolet 
lines, the abundance corrections applicable in the case of 
the infrared OH lines may in practice be smaller than for 
the ultraviolet lines (with the exception that not yet in- 
vestigated non-LTE effects may impact the UV but not 
the IR lines). We note though that no 3D models with 
parameters appropriate for the IR targets at low metallic- 
ities are yet available to verify the continuation of the 3D 
behaviour down to these lower T e g. 

The 3D corrections of abundances derived from [Ol] 
and OH are qualitatively similar in that the corrections 
increase with decreasing metallicity, but the corrections 
for OH are larger than for [Ol]. At [Fe/H] = -2.5, the 
approximate limit of our sample, the OH-based O abun- 
dance is reduced by about 0.5 dex at T c s — 5800 K, and 
by about 0.7 dex at T cB = 6200 K, but the [Oi]-based 
abundances are decreased by only about 0.2 dex. Given 
that we have found OH and [O i] lines to give essentially 
the same abundance with the ID atmospheres, an appli- 
cation of the 3D atmospheres would apparently reopen 
the oxygen problem that we have just closed. A resolu- 
tion of this problem may lie in the replacement of the 
assumption of LTE by non-LTE. The [O i] line is immune 
to non-LTE effects but formation of OH ultraviolet lines in 
non-LTE will lead to weaker lines, and, hence, the need for 
a higher O abundance; Asplund & Garcia Perez estimate 
a rough correction of 0.2 dex but with little [Fe/H] sensi- 
tivity. This correction would not apply to the OH infrared 
lines. An additional factor applicable to both ultraviolet 
and infrared lines is that the OH density may not be at 
its LTE value throughout the atmosphere when granular 
flows are present. 



Since the [O i] and Fe n lines are immune to non-LTE 
effects and the corrections for stellar granulation are rela- 
tive small, we present our abundances with the estimated 
corrections from the 3D model atmospheres as the best in- 
dicator of oxygen abundances for metal-poor dwarfs (Fig. 
|l0| ). After correction for stellar granulation, the index 
[O/Fe] rises with decreasing [Fe/H] to reach the 'low' value 
[O/Fe] ~ 0.3 - 0.4 for [Fe/H] from about -0.7 to -2.0 with 
perhaps an increase to lower [Fe/H]. Certainly, the OH ul- 
traviolet lines as analysed using ID atmospheres indicate 
a rising [O/Fe], but the 3D atmospheres suggest an in- 
creasing correction with decreasing [Fe/H] which cancels 
much of the apparent rise. 

By showing that the OH and [O i] lines give consistent 
abundances when analysed with ID atmospheres, we ap- 
pear to have removed a key part of the oxygen problem. 
Yet, the 3D atmospheres by suggesting lower oxygen abun- 
dances from OH lines than from the [O i] line introduce 
an inverse of the original problem. While consideration 
of non-LTE effects may yet correct this inverse problem, 
the largest potential source of systematic errors in oxygen 
abundance determinations may reside in present realisa- 
tions of stellar granulation. 

5.3. [Ol] and Ol lines 

The permitted triplet of O I lines near 7774 A has been 
analysed extensively in metal-poor dwarfs and giants. 
Analysts of the OH ultraviolet lines have argued that the 
OH and O I lines provide consistent oxygen abundances. 
Other analyses of the O I lines have confirmed that the O I 
lines give 'high' [O/Fe] values for metal-poor stars (e.g., 
Cavallo et al. |1997| ; Mishenina et al. |2000j ). 

As described in Sect. 4.3 we determined oxygen abun- 
dances from a non-LTE analysis of the Ol triplet for 15 
stars. The resulting [O/Fe] values are plotted in Fig. |[ 
As seen the triplet gives [O/Fe] indices similar to the ID 
values from the [O i] line. For five of the 15 stars, we used 
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both the permitted and the forbidden line finding differ- 
ences (permitted - forbidden) ranging from —0.07 to 0.21 
dex and a mean difference of 0.03 dex. This good agree- 
ment exists for analyses using the ID model atmospheres. 



Present indications (Asplund 2001) are that a LTE 
analysis using an equivalent 3D model atmosphere leads to 
a less than 0.1 dex change in the oxygen abundance from 
O I lines. Then, application of 3D model atmospheres re- 
duces the abundance from the [Ol] line more than that 
from the Ol lines. Unless the non-LTE corrections to the 
O I lines are enhanced in the 3D model atmospheres, oxy- 
gen abundances from the O I lines may be ~ 0.2 dex larger 
than from the [O i] line. 

5.4. Oxygen abundances in dwarfs and giants 

Arguments favoring the [O i] - Fe n combination of lines 
as the optimal indicator of the [O/Fe] ratio hold for gi- 
ants as well as dwarfs. Indeed, the run of [O/Fe] versus 
[Fe/H] was considered as well established from [Ol] lines 
in giants well before the appearance of the measurements 
of OH ultraviolet lines in dwarfs. For example, Barbuy 
(|l988p showed that [O/Fe] ~ +0.4 from analyses of metal- 
poor giants with [Fe/H] extending down to —2.5. (Her 
estimates of [Fe/H] are not based on Fell but Fel lines, 
although the few illustrated fits to observed spectra imply 
a good fit to the Sen line near 6300.6 A). Other analyses 
using the [Ol] line gave similar results for [O/F e] (e.g., 
Lambert et al. 1974; Gratton & Ortolani 1986] ; Sneden 
et al. 1991; Kraft et al. 1992). Since most analyses refer- 



enced the stellar oxygen abundance to a solar abundance 
approximately 0.2 dex higher than our adopted value, and 
used a log gf value 0.03 dex smaller than our value, the 
published [O/Fe] should be increased by about 0.2 dex. In 
some cases, a revision should be applied also to account for 
different values of the solar iron abundance. Early analy- 
ses may have adopted a solar abundance loge(Fe) = 7.65 
instead of the now favored value close to loge(Fe) = 7.50. 
Such revisions simultaneously to the solar O and Fe abun- 
dances cancel almost exactly for [O/Fe]. 

Recently, Sneden & Primas ( |2001 ) obtained high- 
quality spectra and analysed the [O i] 6300 A line along- 
with Sen and Fen lines to give [O/Fe] and [O/Sc] with 
[Sc/Fe] ~ 0.0 expected from many studies. After adjust- 
ment for our lower solar oxygen abundance and higher 
g/-value for the forbidden line, their preliminary result 
is [O/Fe] ~ +0.45 with an indication of a slight increase 
with decreasing [Fe/H] to reach [O/Fe] ~ +0.6 at [Fe/H] 
= -2.7 from [O/Fe] ~ +0.45 at [Fe/H] = -1.2. In view of 
the insensitivity of the [O /Fe] results to the adopted model 
atmospheres, these results are unlikely to change appre- 
ciably for alternative choices of effective temperature, and 
surface gravity. 

King ( |2000 ) discussed published analyses of [Ol] and 
Fe II lines in metal-poor giants and obtained almost iden- 
tical results (see his Fig. 5). Much of King's paper is 
concerned with the application of non-LTE corrections 



to published iron abundances derived from Fe I as a way 
to reconcile disagreements over the run of [O/Fe] with 
[Fe/H]. By focussing on [Ol] and Fen lines, non-LTE be- 
comes a moot issue. Insensitivity to non-LTE effects is a 
powerful advantage of using Fe II lines given the prelimi- 
nary nature of extant non-LTE predictions for iron. 

At [O/Fe] ~ +0.5, the giants' [Ol] and Fen lines pro- 
vide results for [O/Fe] comparable to those in Fig. |[ 
The predicted effects of granulation on [O/Fe] for dwarfs 
lower the ratio significantly from values obtained using 
ID model atmospheres. Preliminary corrections give the 
results shown in Fig. [uj|, where [O/Fe] ~ +0.3 for -2.0 < 
[Fe/H] < —1.0 for dwarfs. It is these results that are in 
mild disagreement with those from giants obtained from 
the [Ol] - Fen combination of lines and ID model atmo- 
spheres. To take the final step in the application of the 
[O i] and Fe II lines to a definitive determination of [O /Fe] 
for metal-poor stars, it is necessary to apply 3D model 
atmospheres to analyses of lines in spectra of giants, but 
such 3D model atmospheres do not yet exist. One might 
speculate that effects of stellar granulation on the lines in 
giants will lower [O/Fe] by about 0.2 dex. 

The OH vibration-rotation lines but not the ultravio- 
let lines have been used to ob tain oxygen abundances for 
giants (Melendez et al. 2001). The published results de- 
pend on the choice of atmospheric parameters: [O/Fe] — 
+0.5 or +0.3. These values should be increased by about 
0.13 dex to account for the reduction in the solar oxygen 
abundance. A mean of the two values with the increase of 
0.13 dex gives [O/Fe] ~ 0.5, a result consistent with that 
from the forbidden oxygen line. 

Oxygen triplet lines in giants (logg < 3.0) were ob- 
served and analysed by Cavallo et al. ( 1997| ). Their LTE 
results for [O/Fe] for [Fe/H] < — 1 fall above our points in 
Figjl [O/Fe] ~ +0.8 for [Fe/H] < -1. Mishenina et al. 
( 2000| ) observed a few additional giants and undertook a 
non-LTE analysis for their stars as well as Cavallo et al.'s 
sample. Non-LTE effects as implemented by Mishenina 
et al. reduce [O/Fe] by about 0.1 dex on average, and 
the non-LTE oxygen abundances leave [O /Fe] significantly 
above ours in Fig. |: [O/Fe] ~ 0.37[Fe/H] + 0.11 is a fit 
by eye to Mishenina et al.'s results for giants. Our lower 
solar oxygen abundance implies that the published [O /Fe] 
should be increased by 0.2 dex. A critical analysis of the 
atmospheric parameters for giants including the [Fe/H] in- 
dices might be instructive; we note that for three of the 
four stars in common between Cavallo et al. and Johnson 



& Bolte (2001), an example of a recent spectroscopic anal- 
ysis, there are large differences in the adopted atmospheric 
parameters with the recent choices leading apparently to 
larger [O/Fe] indices. It remains to consider the non-LTE 
effects on the iron lines, or to minimize their influence by 
analysing the Fell lines. In the context of ID model at- 
mospheres, the high oxygen abundance for giants from the 
triplet lines remains the outstanding issue of the oxygen 
problem. 
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6. Concluding remarks 

In this paper, we have demonstrated that an efficient high- 
resolution spectrograph on a very large telescope can now 
provide firm detections of the weak [O i] lines in spectra of 
metal-poor dwarfs. Analysis of the forbidden oxygen lines 
with Fen lines provides an estimate of [O/Fe] which is 
insensitive to the adopted model atmosphere parameters 
T c ff and log g. 

Our results using standard ID model atmospheres 
show [O/Fe] to increase with decreasing [Fe/H] (Fig. ^). 
We have shown that the oxygen problem largely vanishes 
when ID models are applied; in particular we find a sat- 
isfactory agreement between oxygen abundances derived 
from the [Oi] 6300 A line in dwarfs and subgiants and 
abundances obtained from OH ultraviolet and infrared 
lines, and the O I triplet lines. Published results for [O/Fe] 
in metal-poor giants using also the [Oi] - Fen combina- 
tion of lines are compatible with results in Fig. ^| from the 
same combination. Observations of the O I lines in giants 
give LTE abundances and [O/Fe] values larger than our 
[Oi]-based results in Fig. [}. Non-LTE effects reduce the 
discrepancy. 

Application of 3D model atmospheres reintroduces in 
a mild fashion aspects of the oxygen problem. The LTE 
OH-based oxygen abundances for dwarfs appear to be less 
than the [O i]-based results, but the O i-based results are 
appreciably larger than the [O i]-based results in Fig. |l(| 
Non-LTE effects yet to be quantitatively assessed for OH 
and O I lines in the 3D atmospheres may reduce these 
differences relative to the non-LTE insensitive [Oi] line. 
There are also small abundance differences between re- 
sults from 3D atmospheres for dwarfs and ID atmospheres 
for giants. Differences arising from the introduction of the 
3D atmospheres should challenge the confidence of stel- 
lar spectroscopists using the now extensive grids of ID 
model atmospheres. To respond to this challenge theoret- 
ical studies of atmospheric inhomogeneities (granulation) 
must be pursued. There is an urgent need for 3D model 
atmospheres for giants. A parallel effort on observational 
quantification of the effects of granulation on spectral lines 
(and the continuum fluxes) must be carried out: measure- 
ments of line asymmetries and velocity shifts, as well as 
equivalent widths of lines. Lack of understanding of stellar 
granulation may presently be the primary source of sys- 
tematic errors in abundance analyses. Evidence presented 
here indicates that the magnitude of these errors may well 
exceed the errors often confidently attached to abundances 
derived from classical model atmospheres. 

Finally, we note that in describing plots like Fig. [lO] 
there is a tendency to impose a simple linear relation onto 
the points. Instead one might identify a plateau at [O/Fc] 
~ +0.3 for [Fe/H] from about -0.7 to -2.0 with a rise to 
higher [O/Fe] at lower [Fe/H]. Furthermore, not all stars 
may lie on a single relation. Although much of the ap- 
parent scatter about a mean relation is attributable to 
measurement errors, there are a few stars that most prob- 
ably fall off a mean relation. HD 103273 with a lower than 



average oxygen abundance is a halo star known to be un- 
usually deficient in Mg, Si, and Ca (a-elements) (Nissen & 
Schuster 1997), and a slight oxygen deficiency fits this pat- 
tern. G 75-31 may be another example of an cv-poor star. 
HD 78558 is a thick-disk star with an above average [a/Fe] 
and, hence, its unusually high [O/Fe] is not unexpected. 
The star has a small mean galactocentric distance accord- 



ing to Nissen & Edvardsson (1992). In order to investigate 
the reality of this possible cosmic scatter of [O/Fe] at a 
given value of [Fe/H] a much larger sample spanning stars 
with different Galactic orbits should be studied. 
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